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DYNAMIC RESPONSE OF A FAMILY OF AXISYMMETRIC HAMMERHEAD 


MODELS TO UNSTEADY AERODYNAMIC LOADING 

By Robert C. Robinson, Phillip R. Wilcox, 
Bruno J. Gambucci, and Robert E. George 

Ames Research Center 


SUMMARY 


The effects of boattail angle and diameter ratio on the unsteady aerody- 
namic loading on hammerhead launch vehicles were studied in wind-tunnel tests 
of 13 related dynamic models using the partial -mode model technique. The test 
Mach numbers ranged from 0.80 to 2.50 and the Reynolds numbers from 3-2X10 6 to 
4.5X10 6 based on maximum diameter of the models. It was found that the 
unsteady aerodynamic loading was the result of two phenomena: buffeting due 

to separated flow and dynamic instability due to fluctuations between 
separated and attached flow. 


INTRODUCTION 


Many space payloads have a diameter greater than that of the final rocket 
stages used in launching, thus the launch configurations are "hammerhead" 
shapes. The geometrical transition between the payload and the rocket may 
cause a region of separated flow and, consequently, large fluctuations in the 
local pressure, such as those discussed in references 1 and 2. If those pres- 
sure fluctuations were correlated over a large area, a coupling between pres- 
sure and motion could develop that would transfer aerodynamic energy to the 
bending modes or rigid body dynamics of the complete vehicle. Fluctuations 
between separated and attached flow due to changes in angle of attack can, 
with the proper phasing, also transfer energy to the vehicle. Such phenomena 
have been observed in wind-tunnel tests of dynamic models, and the measured 
response and aerodynamic damping of several models have been reported in ref- 
erence 3 . Data from two of those configurations are included in this report. 
In reference 4 it is shown that the response of a free-oscillation model can 
be well predicted from pressure fluctuations measured on a static model if 
aerodynamic damping is included in the calculations. 

The purpose of the present series of tests was to investigate the effect 
of a systematic variation of diameter ratio and boattail angle on the dynamic 
stability of partial mode models of hammerhead vehicles. Another objective 
was to correlate the stability of the models with the stability criterion of 
reference 5 - 
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half -power bandwidth, Hz 

rms resp onse moment at the node 
qSD 

maximum diameter 
base diameter 

model resonant frequency, Hz 
free -stream Mach number 
free -stream dynamic pressure 
maximum cross sectional area 
angle of attack 
boattail angle 

ratio of damping to critical damping 
tare damping ratio 


amplitude spectral density, degrees/s/Hz 


amplitude spectral density at f = 5 Hz 


MODELS AND TEST EQUIPMENT 


Models 

In the partial -mode model technique used in the present tests, rigid mod- 
els are designed to simulate the portion of the full-scale vehicle where the 
important aerodynamic forces are expected to occur, as explained in some 
detail in references 4 and 6. The aerodynamic forces that contribute most to 
buffeting and instability of hammerhead shapes have been found to occur in the 
boattail area. Models for these tests were designed to simulate a mode in 
which the most forward node is behind the boattail. The portion of the model 
behind the node is attached rigidly to the support and the response of the 
active part of the model is measured at the node. 

Two basic models with the same nose shape were used. One model had a 
maximum diameter of 11.19 inches (28. 42 cm) and a base diameter of 7*00 inches 
(l7*78 cm). The other model had a maximum diameter of 12.00 inches (30.48 cm) 
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and a base diameter of 6.00 inches (15-24 cm). The profile and dimensions of 
the larger model are shown in figure l(a); the several combinations of 
boattail angle and diameter ratio that were tested are shown in figure l(b). 

Boattail angle and diameter ratio were changed by the use of "slip-ons" 
which were slid over the model base and brought to the correct position on the 
model where they were glued in place. Figure 2 shows the two types of slip- 
ons. With this technique changing, the model configuration was rapid and 
inexpens ive . 

The construction of both models was quite similar except that the rear 
cylindrical section of the large model, with a diameter ratio of 2 , was metal 
because of the diameter limitations of the balance. Magnesium was selected 
for the cylinder to keep the model weight to a minimum. In general, the mod- 
els were of lightweight sandwich -type construction. The inner shell was made 
of four layers of 0.004 inch ( 0.1016 mm) fiber glass cloth, and the outer 
shell of two layers of the same material; the filler was polyurethane foam 
having a specific weight of 2 lb/cu ft (31-4 N/m 3 ) . Polyester resin was used 
sparingly as the binder. The smaller model, with a diameter ratio 1 . 6 , was 
made entirely of the sandwich material except for a thin, short aluminum 
cylinder bonded to the aft end of the inner shell for mounting the model on 
the balance. 

A 1-inch strip of number 40 carborundum grit, for inducing boundary-layer 
transition, appears as a dark band on the nose of the model in figure 3« The 
choice of grit size was based on the investigations reported in references 7 
and 8 . 


Wind Tunnels 

The tests were conducted in the Ames l4-Foot (4.27 meter) Transonic Wind 
Tunnel, 11-Foot (3-35 meter) Transonic Wind Tunnel, and 9“ hy 7“Foot 
(2.74 X 2.13 meter) Supersonic Wind Tunnel. The Mach number was varied from 
0.80 to 1.40 in the transonic wind tunnels and from 1-55 to 2.50 in the super- 
sonic wind tunnel. The corresponding Reynolds number range was 3.5 X 10 6 to 
4.5X10 6 based on the maximum diameter of the models. 

In all three facilities, the models were mounted on a sting attached to a 
movable strut. In the transonic wind tunnels, both the strut and the model 
plane of oscillation were vertical (figs. 3 and 4); whereas in the supersonic 
tunnel, although the plane of oscillation was vertical, the strut operated in 
a horizontal plane (fig. 5 )* 


Balance 

Details of the free -oscillation balance used In the tests are shown in 
figure 6 (b). The model is secured to the model mount which is attached to a 
short cantilever resonant spring. The resonant spring serves a dual purpose: 
It fixes the model axis of rotation at the desired node line, and the output 
from Its strain gages is the time history of the model motion about the center 
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of rotation. The air brake is used to limit the model deflection and thus 
keep the resonant spring bending stresses at a safe level. The maximum deflec- 
tion of the model is controlled by a limit switch that energizes an air valve 
to operate the air brake. Figure 6(a) shows the relation of the fixed and 
active parts of a typical model to the sting and balance. 


Ins trument at i on 

The data from these tests can be classified into two groups: wind-off 

data used in evaluating the dynamic characteristics of the models and support 
systems, and wind-on data which were analyzed to evaluate the effects of con- 
figuration changes on the aerodynamic forces. Figure 7 shows the instrumenta- 
tion used to record the model response during both wind-off and wind-on tests. 
In the "operate" position the standard termination is AC coupled to eliminate 
any bias due to static aerodynamic loads, thus making the entire dynamic range 
of the recording system usable for the dynamic response of the models. In the 
"calibrate" position DC is passed, permitting static calibration of the strain 
gage. The precise AC calibration signal, which was recorded simultaneously on 
al 1 channels, was used as a voltage and frequency reference in the analysis of 
the recorded data. 

Additional instrumentation required in the wind-off tests included a fre- 
quency response analyzer and a system for exciting the model. Figure 8 shows 
the arrangement for measuring the frequency response of a model. The output 
of the force gage in the impedance head and the output of the strain gage in 
the balance were sent to the frequency response analyzer which produced a log 
plot of amplitude ratio versus frequency. The same system was used in measur- 
ing the sting characteristics except that both inputs to the analyzer were 
taken from the impedance head. The dynamic characteristics of each model were 
checked before and after each wind-on test by recording a free oscillation 
decay. In the earlier tests the models were excited by hand tapping, but an 
improved technique was developed (fig. 9) in which a step input was produced 
when a weight hung from an elastic strap was released by cutting the connect- 
ing wire. This provided both a well defined input to the model and a check on 
the static calibration of the strain gage. 


ESTIMATION OF TAKE DAMPING 


To evaluate the effects of configuration changes and test conditions on 
the aerodynamic damping, it is essential to have a good approximation to 
single -degree -of -freedom motion and a nearly constant, low level of tare damp- 
ing. If the model frequency is too near a support frequency, the model 
response will contain two frequencies, and the effective damping of the model 
will be changed. To avoid this condition either the model frequency or the 
support frequencies may be altered. The model frequencies were changed by 
adding mass or by changing to a different resonant spring. The support char- 
acteristics were changed by tethering the sting (as shown in fig. 4), by 
suspending a weight from the sting, or by changing to a sting of different 
length. 
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A short sting with a high first -bending frequency was used for most of 
the tests in the 11-foot transonic wind tunnel. The frequency response for 
this sting (fig. 10) shows a first -bending frequency of 38 Hz, well above the 
desired model frequency of 20 Hz (fall-scale frequency of 2.0 Hz). A fre- 
quency response plot for one of the larger models mounted on this sting is 
shown in figure 11 where it can be seen that at the model frequency of 19-2 Hz 
a single degree of freedom was very well approximated. At the sting frequency 
the deviation from the single -degree -of -freedom curve is only about 0.5 per- 
cent of the model’s peak response. The corresponding free -oscillation decay 
for this model (fig. 12(a)) is smooth with f n = 20 Hz and £ =• 0 . 016 . However, 
when the smaller models were mounted on the short sting, the small increase in 
model resonant frequency resulted in an uneven decay, as shown in figure 12 (b). 
Analysis of this decay by the technique reported in reference 9 showed the 
model frequency to be 21 Hz with £ = O.O 35 . The interfering frequency was 
found to be 23.5 Hz. This unknown support mode is barely discernible in the 
sting frequency response (fig. 10 ) . Fortunately the static instability of 
these configurations about the moment center lowered the wind-on resonant fre- 
quency below 20 Hz. Consequently interaction with the support was small as 
shown by the typical amplitude spectra of the wind-on response presented in 
figure 13 . Because of the decreased frequency with wind on, it is estimated 
that the wind-on tare damping of the models which had wind-off frequencies 
above 20 Hz was about the same as for models with wind-off frequencies below 
20 Hz, namely 0.012 ^ t, £ 0.015. Exceptions to this are models 3 and H which 
were tested on a different sting and had a very low level of tare damping. 

Table I lists the wind-off frequency and damping and the average wind-on 
frequency for all the models . 


AERODYNAMIC EFFECTS ON MODEL RESPONSE 


In the wind-on tests, the data acquisition procedure consisted in record- 
ing a 60-second time history of the model at fixed Mach numbers and angles of 
attack. Changes in Mach number and attitude were made with the models free to 
oscillate about the node line so that any large response confined to a narrow 
range of M or a could be detected. The model response was monitored contin- 
uously on an oscilloscope, and if any significant increases in amplitude were 
noted, data were recorded at those points in addition to those scheduled. If 
the model response exceeded a preset level, the brake was automatically actu- 
ated. In some cases response of the model was great enough to cause brake 
actuation before a record long enough for analysis could be obtained. 

The random character of the model response is illustrated by the two typ- 
ical time histories in figure l4. The amplitude distribution of the absolute 
values of the responses was determined by the technique described in refer- 
ence 10. The data in figure 15 follow the normal distribution and therefore 
the amplitude distribution can be completely characterized by the root mean 
square. The rms of the strain-gage output was obtained by integrating over 
the full length of each record and was then reduced to a response -moment 
coefficient, C mrms - Figures l 6 through 19 show the variation of Cmrms ‘with 
Mach number at a = 0° and 2 ° for all the models. 


5 



II 


Six points for which the response was exceptionally large have been 
labeled "unstable" on figures 18 and 19 . However, dynamic stability cannot be 
judged on the basis of amplitude alone as a large response can result from 
either negative aerodynamic damping or increased buffeting. A measure of the 
aerodynamic damping can be obtained from amplitude spectra such as those shown 
in figure 13 . It can be shown (ref. 11 ) that for small damping the half -power 
bandwidth of the spectrum peak is related to the damping ratio by the follow- 
ing equation: 


B - 2 £f n 

A "bandwidth, measurement which shows the net damping (tare damping plus 
aerodynamic damping) to he less than the tare damping is an indication of 
negative aerodynamic damping. 

The variation of bandwidth with Mach number for the D/d = 1.6, S — 17° 
model at a = 0° and 2° is shown in figure 20. The tare damping converted to 
bandwidth is also shown for comparison. At zero angle of attack, the net 
damping is much less than the estimated tare damping at M = 1.00, confirming 
t ha t the large response for this test condition is due to dynamic instability. 
The bandwidth at a = 2° shows a sharp drop at M = O .98 and O. 99 , indicating 
that a very narrow range of instability may exist at a = 2 ° and M O. 985 . 

The D/d = 2.0, S = 20° model was definitely unstable at a = 2° as well 
as at cl - 0 ° (fig. 2l) . At a = 2 ° and M = l.l 8 the total damping approaches 
zero, whereas at the other Mach numbers it is several times greater than the 
tare damping for this angle of attack. For a = 0°, the response was so vio- 
lent at M = 1.25 that a record long enough for analysis could not be obtained; 
however, the trend of the bandwidth data is toward -zero damping for this test 
condition. The variation of bandwidth with angle of attack in figure 22 shows 
that instability occurred also at M = 1.20 and a = 1.4°. 

Although the D/d =1.6, 5 = 20° model was unstable at M = 1.20, a = 0° 
in the l4-foot wind-tunnel tests, as reported in reference 1 , it was stable in 
the 11-foot wind-tunnel tests. Figure 18(d) shows that the response at 
M = 1.20 was very low and bandwidth measurements indicated that the 
aerodynamic damping was positive for this model at all the test conditions. 

To show more clearly the interrelated effects of diameter ratio, boattail 
angle, and Mach number, isometric plots were made showing Cmpms as a : ^ unc " 
tion of those parameters. Figure 23 shows such plots for diameter ratios of 
1 . 25 , 1.60, and 2.00 at a = 0°. For D/d = 1.25 the maximum response occurs 
between M = 0.90 and M = 1.00 forming a broad level ridge. The effect of 
boattail angle is negligible. At D/d = 1.60 (fig. 23 (b)) there are two 
ridges. One is relatively low and broad with its maximum at M = 0.95 for all 
values of 5. It is considered to be the result of a high level of buffeting 
due to separated flow over the boattail. The other ridge, narrow and of large 
amplitude, is formed by the three unstable points. The two ridges coincide at 
5 = 10° and M = O. 95 . As S is increased to 17°, the Mach number for insta- 
bility increases to 1.0. There is a buffeting ridge between M = 0.90 and 
M = 1.00 for D/d = 2.00 also (fig. 23 (c)), and one point of instability is at 
M = I. 25 . The larger Mach number separation between the buffeting maximum and 
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the unstable point for this diameter ratio makes it more apparent that two 
different phenomena are involved. One is the buffeting due to separated flow 
and the other is the instability due to fluctuations between separated and 
attached flow. A similar plot for the D/d = 2.0 models at a - 2° would show 
a smaller Mach number separation between the unstable point and the buffeting 
ridge . 

A comparison of the three parts of figure 23 shows the differing effects 
of changes in boattail angle and diameter ratio. It is evident that although 
the boattail angle is a determining factor in the dynamic stability of a con- 
figuration, it has only a small effect on the buffeting level. Diameter ratio 
however, has an important influence on both dynamic stability and buffeting 
level. There appears to be a minimum D/d below which instability does not 
occur for a given boattail angle, and there is a significant increase in 
buffeting with increasing diameter ratio. 

The stability criterion of reference 5 is based on the pressure gradient 
in the boattail region computed from a theoretical pressure distribution and 
classifies configurations having maximum adverse pressure gradients between 

0.2 and 1.8 as being "separated-unstable. " The maximum adverse pressure gradi 
ents for the present 13 configurations varied from O .33 to 2 . 23 . The pressure 
gradients for the four dynamically unstable configurations fall in the 
separated-unstable range; however, eight other configurations in that range 
were stable, on the basis of data from these tests. 


CONCLUSIONS 


Results from wind-tunnel tests of partial -mode models of 13 related 
hammerhead conf igurations lead to the following conclusions: 

1. Unsteady aerodynamic loading on launch vehicles with hammerhead 
payloads can result from two phenomena: buffeting caused by separated flow 

and dynamic instability due to fluctuations between separated and attached 
flow. 


2. The Mach number and angle of attack at which dynamic instability 
occurs are influenced by both boattail angle and diameter ratio. 

3 . The buffeting level increases with diameter ratio but is insensitive 
to changes in boattail angle. 

k. The stability criterion of NASA SP 8001 is conservative for 
configurations similar to those studied in these tests. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 9^035, Dec. 28, 1967 
124-08-04-09-00-21 
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TABLE I. - FREQUENCY AMD DAMPING CHARACTERISTICS OF THE MODELS 


Model 

D/d 

5 

Wind 

Wind- 

off 

Wind-on 

tunnel 

f n > Hz 

5 

Average f n , Hz 

i 

2.0 

30 ° 

11-foot 
9- by 7 -foot 

20.0 

17-9 

0.016 

0.009 

16.5 

16.0 

2 

2.0 

25 ° 

11 -foot 
9 - "by 7 -foot 

19-5 

17.8 

0.012 

0.007 

16.1 

15-7 

3 

2.0 

20° 

11-foot 

19.4 

0.002 

16.4 

4 

1.6 

30 ° 

11-foot 
9- "by 7-foot 

19-2 

17.8 

0.015 

0.009 

16.5 

15.6 

5 

1*6 

25 ° 

11-foot 

21.0 

0.035 

19.2 

6 

(ref. 3) 

1.6 

20° 

11 -foot 
l 4 -foot 

21.0 

23.3 

0.018 

0.017 

18.7 

7 

1.6 

17 ° 

11-foot 

21.0 

0.034 

18.5 

8 

1.6 

15 ° 

l 4 -foot 

19.7 

0.019 

— 

9 

(ref. 3) 

1.6 

10° 

l 4 -foot 

24.6 

0.018 

— 

10 

1.25 

30 ° 

11 -foot 

I8.7 

0.014 

16.5 

9“ "by 7 -foot 

17.6 

0.014 

15.5 

11 

1.25 

20° 

11-foot 

18.7 

0.004 

16.7 

12 

1.25 

17 ° 

11-foot 

20.5 

0.018 

18.5 

13 

1.10 

30 ° 

11 -foot 

I8.7 

0.014 

16.2 

9- "by 7-foot 

17.3 

0.007 

15.0 
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(a) Slip-on for changing model boattail angle. 







Figure 3 .- A model mounted in the 14-foot transonic wind tunnel 
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Figure 4 . - A model mounted in the Ames 11 -Foot Transonic Wind Tunnel. 
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Figure 


(b) Details of the free-oscillation balance. 

6.- The method of mounting the models on the sting. 















A-35375. 1 

Figure 8 .- Shaker installation for wind -off tests. 
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Figure 9-- The technique used to produce a step input for initiating free 

oscillations . 
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Figure 1 6 .- Variation of the moment coefficient with Mach number for the D/d = 1.10, 8 = 30° model. 









Figure l8.- Variation of moment coefficient with Mach number for models with 

D/d = 1.6. 
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(a) 5 = 20° 

Figure 19*" Variation of moment coefficient with ]VIach number for the models 

with D/d = 2.0. 
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(b) D/d = 1.6 


Figure 23 .- Continued. 
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